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ABSTRACT 

We analyse the spatial clustering properties of the RASSl Bright Sample, 
an X-ray flux-limited catalogue of galaxy clusters selected from the southern 
part of the ROSAT All-Sky Survey. The two-point correlation function ^(r) 
of the whole sample is well fitted (in an Einstein-de Sitter model) by the 
power-law ^ = (r/ro)~'^, with vq = 21.5tli Mpc and 7 = S.ll^gj^ 
(95.4 per cent confidence level with one fitting parameter). We use the RASSl 
Bright Sample as a first application of a theoretical model which aims at 
predicting the clustering properties of X-ray clusters in flux- limited surveys for 
different cosmological scenarios. The model uses the theoretical and empirical 
relations between mass, temperature and X-ray cluster luminosity, and fully 
accounts for the redshift evolution of the underlying dark matter clustering 
and cluster bias factor. The comparison between observational results and 
theoretical predictions shows that the Einstein-de Sitter models display too low 
a correlation length, while models with a matter density parameter J7o,„ = 0.3 
(with or without a cosmological constant) are successful in reproducing the 
observed clustering. The dependence of the correlation length tq on the X- 
ray limiting fiux and luminosity of the sample is generally consistent with the 
predictions of all our models. Quantitative agreement is however only reached 
for rJom = 0.3 models. The model presented here can be reliably applied to 
future deeper X-ray cluster surveys: the study of their clustering properties 
will provide a useful complementary tool to the traditional cluster abundance 
analyses to constrain the cosmological parameters. 

Key words: cosmology: theory ~ galaxies: clusters - large-scale structure of 
Universe - X-rays: galaxies 



1 INTRODUCTION 

Galaxy clusters play an important role in the models for structure formation based on the gravitational instability 
hypothesis. They are the most extended gravitationally bound systems in the Universe. For this reason the study of 
their properties is a useful tool to constrain the parameters entering in the definition of the cosmological scenarios. In 
particular, their abundance and spatial distribution (also as a function of redshift) have been used to obtain estimates 
of the mass fluctuation amplitude and of the density parameter f2o (e.g. Eke, Cole & Frenk 1996; Viana & Liddle 
1996; Mo, Jing & White 1996; Oukbir, Bartlett & Blanchard 1997; Eke et al. 1998; Sadat, Blanchard & Oubkir 1998; 
Viana & Liddle 1999; Borgani, Plionis & Kolokotronis 1999; Borgani et al. 1999). In the past years many different 
groups have compiled deep cluster surveys in the optical band, which have been used to compute the clustering 
properties of galaxy clusters. The flrst results showed that clusters are strongly correlated, with a correlation length 
ro ~ 20 — 25 Mpc, a factor 4-5 larger than that obtained for local galaxies (e.g. Bahcall & Soneira 1983; Postman, 
Huchra & Geller 1992). Here h represents the Hubble constant Ho in units of 100 km s"'^ Mpc"'^. Sutherland (1988) 
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suggested the existence of a possible strong effect due to the spurious presence of galaxies acting as interlopers (see 
also Dekel et al. 1989; van Haarlem, Frenk & White 1997). New analyses of optical catalogues, taking into account 
this projection effect (Dalton et al. 1992; Nichol et al. 1992; Dalton et al. 1994; Croft et al. 1997), led to a smaller 
value for the cluster correlation length (ro « 13 — 18 Mpc). 

A way to overcome the projection problem is the use of data obtained in the X-ray region of the spectrum. In 
fact, in this band, galaxy clusters have a strong emission produced by thermal bremsstrahlung, which allows to detect 
them also at high redshifts. Starting from the eighties, different space missions produced extended cluster catalogues 
which have been mainly used to compute their X-ray luminosity function. In particular, the ROSAT satellite provided 
a good opportunity to build a reliable all-sky survey, which was performed in the soft (0.1 - 2.4 keV) X-ray band. 
First studies of the clustering properties in small samples of X-ray selected galaxy clusters have been performed by 
Lahav et al. (1989), Nichol, Briel & Henry (1994) and Romer et al. (1994). More recently the ROSAT data have 
been correlated with the Abell-ACO cluster catalogue (Abell, Corwin & Olowin 1989) to produce the X-ray Brightest 
Abell Cluster sample (XBACs; Ebeling et al. 1996), for which estimates of the two-point correlation function have 
been recently obtained (Abadi, Lambas & Muriel 1998; Borgani, Plionis & Kolokotronis 1999). The corresponding 
values for ro are in the range ro « 20 — 26 Mpc. A smaller amplitude of the correlation function is obtained from 
the preliminary analyses of the REFLEX sample (Collins et al. 1999), which is also obtained by the ROSAT All-Sky 
Survey (RASS) data. 

In this paper we estimate the clustering properties for the RASSl Bright Sample (De Grandi et al. 1999), which 
is another X-ray cluster catalogue obtained from the RASS. In this case the clusters are spectroscopically searched 
in a preliminary list of candidates produced by correlating the X-ray data with regions of galaxy overdensity in the 
southern sky. In this way, the resulting catalogue is not affected by the selection biases present in the Abell-ACO 
cluster catalogue. The RASSl Bright Sample is used to test a theoretical model for the correlation function of X-ray 
clusters in Hux- limited samples (see also Moscardini et al. 1999). This model makes use of the technique introduced 
by Matarrese et al. (1997) and Moscardini et al. (1998), which allows a detailed modelling of the redshift evolution of 
clustering, accounting both for the non-linear dynamics of the dark matter distribution and for the redshift evolution 
of the bias factor. A characteristic feature of this technique is that it takes into full account light-cone effects, which 
are relevant in analysing the clustering of even moderate redshift objects (see also Matsubara, Suto & Szapudi 1997; 
de Laix & Starkman 1998; Yamamoto & Suto 1999). 

The plan of the paper is as follows. In Section 2 we summarize the characteristics of the RASSl Bright Sample 
used in the following clustering analysis. In Section 3 we discuss the method used to compute the observational 
two-point correlation function for the RASSl Bright Sample and we present the results. In Section 4 we introduce 
our theoretical model to estimate the correlations of the X-ray clusters in the framework of different cosmological 
models and we compare our predictions to the observational results. Conclusions are drawn in Section 5. 



2 THE SAMPLE 

The RASSl Bright Sample (De Grandi et al. 1999), contains 130 clusters of galaxies selected from the first processing 
of the ROSAT All-Sky Survey (RASS) data (Voges 1992). This sample was constructed as part of an ESQ Key 
Programme (Guzzo et al. 1995) aimed at surveying all southern RASS candidates, which is now known as the 
REFLEX cluster survey (Bohringer et al. 1998; Guzzo et al. 1999). The identification of RASS cluster candidates was 
performed by means of different optically and X-ray based methods. First, candidates were found as overdensities 
in the galaxy density distribution at the position of the X-ray sources using the COSMOS optical object catalogue 
(e.g. Heydon-Dumbleton, Collins & MacGillivray 1989). Then, correlating all RASS sources with the AGO cluster 
catalogue, and, finally, selecting all RASS X-ray extended sources. X-ray fluxes were remeasured using the steepness 
ratio technique (De Grandi et al. 1997) , speciflcally developed for estimating fluxes from both extended and pointlike 
objects. A number of selections aimed at improving the completeness of the final sample lead to the RASSl Bright 
Sample. Considering the intrinsic biases and incompletenesses introduced by the X-ray flux selection and source 
identiflcation processes, the overall completeness of the sample is estimated to be ^ 90 per cent. The RASSl Bright 
Sample is count-rate-limited in the ROSAT hard band (0.5 - 2.0 keV), so that due to the distribution of Galactic 
absorption its effective flux limit varies between 3.05 and 4 x 10~^^ erg cm~^ s~^ over the selected area. This covers 
a region of approximately 2.5 sr within the Southern Galactic Cap, i.e. 5 < 2.5° and b < —20°, with the exclusion of 
patches with RASS exposure times lower than 150 s and of the Magellanic Clouds area. The exact sky map covered by 
the sample is shown in Figure 2 of De Grandi et al. (1999). The redshift distribution for our whole sample is presented 
in the left panel of Figure ^ while the X-ray luminosity Lx as a function of the redshift for each cluster is shown in 
the right one. It is possible to notice that 66 per cent of the clusters have z; < 0.1 but the redshift distribution has a 
tail up to 2: ~ 0.3. 
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Figure 1. Left panel: the redshift distribution of the RASSl Bright Sample. Right panel: the X-ray luminosity Lx (in the 
ROSAT hard band 0.5 - 2.0 keV) vs. the redshift z for each cluster of the sample. The solid line shows the X-ray luminosity 
corresponding to the limiting flux Siim = 3 X 10~^^ erg s~^ cm"'^ in an Einstein-de Sitter model. 



3 THE 2-POINT CORRELATION FUNCTION 

Before computing the clustering properties of our sample, we have to derive for each cluster the comoving radial 
distance r from the observer, given the redshift of each source. To this goal we use the standard relation (neglecting 
the effect of peculiar motions) 

r{z)^ J S ( V|»o7z| / dz'Ul + z'Y {l + Q,o,,J)-z' {2 + z')Q,oa\ ), (1) 

where Qoiz = 1 — fiom — ^oa, with S7om and fioA the density parameters for the non-relativistic matter and cosmological 
constant components, respectively. In this formula, for an open universe model, flon > and S{x) = sinh(a;), for a 
closed universe, Qon < and S{x) = sin(a;), while in the Einstein-de Sitter (EdS) case, flon = and S(x) = x. 

To compute the spatial two-point correlation function ^(r) we adopt both the Landy & Szalay (1993) estimator, 

_ NrjNr - 1) DDjr) - 1 DR{r) 

' N,{N, - 1) RR{r) i?i?(r) ^ ' ^' 

and the Davis & Peebles (1983) estimator, 

^^'^^ " V, - 1 i?7?(r) ^•^^ 

In the previous formulas TV^ is the number of random points and Nc that of clusters, DD is the number of distinct 
cluster-cluster pairs, DR is the number of cluster-random pairs and RR refers to random-random pairs with separation 
between r and r + Ar. The random catalogue contains a number of sources 1,000 times larger than the real catalogue 
(i.e. Nr = lOOOA^c). To generate this sample we have extracted randomly coordinates from the surveyed area (see 
Figure 2 in De Grandi et al. 1999), assigning to each position a random flux drawn from the observed number 
counts (Figure 8 in De Grandi et al. 1999). We decided to retain the source in the catalogue if its flux is larger 
than the flux limit at the choosen position. We adopt two different methods to assign the random redshifts: in the 
first we scramble the observed redshifts of the clusters in the sample; in the second we generate them randomly 
from the observed redshift distribution binned in intervals of 0.01 in z. The results obtained with these two different 
methods are practically indistinguishable. The same happens also if we use the two previous estimators for the 
two-point correlation function (eqs.Pj^. The errorbars have been estimated by using the bootstrap method with 50 
resamplings. We find that the errors obtained in this way are in many cases larger than \/3 times the Poissonian 
estimates, which are often used as an analytical approximation of the bootstrap errors (Mo, Jing & Borner 1992). 
This is particularly true at small separations. 
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Figure 2. Left panel: the two-point correlation function of the whole RASSl Bright Sample. Filled circles, open circles and 
open squares refer to results obtained by assuming (f2omif2oA) = (liO), (Slomif^oA) = (0.3,0) and (nomiS^OA) = (0.3,0.7), 
respectively. The corresponding best-fit relations are shown by heavy solid line, light solid line and dotted line, respectively. 
The 1-(T errorbars obtained by the bootstrap method are shown only for the EdS case for clarity; their sizes in the other cases 
are similar. Right panel: confidence contours (68.3, 95.4 and 99.73 per cent confidence levels) of the fitting parameters ro and 7 
for the two-point correlation function of the whole RASSl Bright Sample. Heavy solid lines and filled circle refer to confidence 
levels and to the best fitting value obtained for the EdS case; light solid lines and open circle show the corresponding results 
for (Qomjf^oA) = (0.3,0) while dotted lines and open square refer to (r2om, ^^oa) = (0.3,0.7) model. 



In the left panel of Figure ^ we show the correlation function computed for the whole catalogue. We present 
results obtained by using both an EdS model and two models with Slom = 0.3 (with and without cosmological 
constant). The results are quite similar and the differences are always small. 

The correlation function has been fitted by adopting the power-law relation 

i{r) = (r/ro)-^ . (4) 

The best-fit parameters have been obtained by using a maximum likelihood estimator based on Poisson statistics and 
unbinned data (Croft et al. 1997; see also Borgani, Plionis & Kolokotronis 1999). Unlike the usual ^'^-minimization, 
this method allows to avoid the uncertainties due to the binsize, to the position of the bin centres and to the bin 
scale (linear or logarithmic). 

To build the estimator, it is necessary to estimate the predicted probability distribution of cluster pairs, given 
a choice for the correlation length ro and the slope 7. By using all the distances between the cluster-random pairs, 
we can compute the number of pairs g{r)dr in arbitrarily small bins dr and use it to predict the mean number of 
cluster-cluster pairs h{r)dr in that interval as 

Kr)dr^^^^[l + i{r)]g{r)dT , (5) 

where the correlation function ^ is modelled with a power-law as in eq.(^). Actually the previous equation holds only 
for the Davis & Peebles (1983) estimator [eq.(^] but, since we obtain very similar results using different estimators, 
we can safely apply it here. Now it is possible to use all the distances between the A'^p cluster-cluster pairs data to 
build a likelihood. In particular, the likelihood function C is defined as the product of the probabilities of having 
exactly one pair at each of the intervals dr occupied by the cluster-cluster pairs data and the probability of having 
no pairs in all the other intervals. Assuming a Poisson distribution, one finds 

£ = cxp[— /i(r)dr]/i(r)dr JJ^ exp[— /i(r)dr] , (6) 

where j runs over all the intervals dr where there are no pairs. It is convenient to define the usual quantity S — -~2\nC 
which can be written, once we retain only the terms depending on the model parameters ro and 7, as 
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5=2 / h{r)dr ~ 2^\n[h{ri)] . (7) 

The integral in the previous equation is computed over the range of scales where the fit is made. We will adopt 5 and 
80 Mpc for rmin and rmax, respectively. 

By minimizing 5* one can obtain the best-fitting parameters ro and 7; the confidence levels are defined by 
computing the increase As with respect the minimum value of S and assuming a distribution for As. 

By applying this maximum likelihood method to the RASSl Bright Sample with the assumption of an EdS model, 
we find ro = 21.514 4 ti~^ Mpc and 7 = 2.1llQ ge (95.4 per cent confidence level with one fitting parameter). Since the 
redshift distribution is shallow, the values obtained in other cosmologies are quite similar: for (f7om,fioA) = (0.3,0) 
we find ro = 21.4+^;g /i"^ Mpc and 7 = 2.17^111, while for (riom,noA) = (0.3,0.7) we find ro = 22.1+|^;7 h"^ Mpc 
and 7 = 2.061'lQ gg. The best-fit relations are also shown in Figure ^ Notice that a x^-minhnization procedure gives 
similar results, but with larger errorbars. 

In the right panel of the same figure we show the contour levels corresponding to As equal to 2.30, 6.31 and 11.8. 
Assuming that As is distributed as a distribution with two degrees of freedom, they correspond to 68.3, 95.4 and 
99.73 per cent confidence levels, respectively. Notice that by assuming a Poisson distribution the method considers 
all pairs as independent, neglecting their clustering. Consequently the resulting errobars can be underestimated (see 
the discussion in Croft et al. 1997). 

Our results are somewhat larger than those derived by Romer et al. (1994) who found ro = 13 — 15 h"^ Mpc 
by analysing a sample of galaxy clusters also selected from the RASS in a similar region of sky (22'' < RA < S*", 
—50° < 5 < 2°, |b| > —40°). A partial explanation of this difference is related to the deeper limiting fiux (Sum — 10~^^ 
erg s^^ cm~^ in the 0.1 - 2.4 keV band) of their catalogue: as we will discuss in the next section, the correlation 
length is expected to depend on the characteristics defining the surveys, such as their limits in fiux and/or luminosity. 
However we have to remind that this early sample was derived drawing on X-ray information from the ROSAT 
standard analysis software (SASS), which was not optimazed for the analysis of extended sources (for a more detailed 
discussion see e.g. De Grandi et al. 1997), and this source of incompleteness was not included in the analysis of Romer 
et al. (1994). Moreover in their analysis the sample sky coverage (i.e., the surveyed area as a function of the fiux 
limit) was not discussed. 

Previous analyses of the XBACs sample, which is a flux-limited catalogue of X-ray Abell clusters with a limiting 
flux Slim = 5 X 10~^^ erg s~^ cm~^ in the 0.1 - 2.4 keV band, gave compatible amplitudes for the correlation function: 
ro = 21.ll2.3 Mpc (Abadi, Lambas & Muriel 1998) and ro = 26.0l4;J, 'i"^ (Borgani, Plionis & Kolokotronis 
1999; errorbars in this case are 2-a uncertainties). Preliminary analyses of the clustering properties of the REFLEX 
sample (Collins et al. 1999; Guzzo et al. 1999), which has a limiting flux Sum = 3 x 10~^^ erg s~^ cm~'^ in the 0.1 - 
2.4 keV band, lead to a smaller correlation length (ro — 18 Mpc). Also in this case, the discrepancy is probably 
a consequence of the deeper limiting flux. 



4 COMPARISON WITH THEORETICAL MODELS 
4.1 Structure formation models 

In the following analysis we consider flve models, all normalized to reproduce the local cluster abundance. In particular 
we will adopt the normalizations obtained by Eke, Cole & Frenk (1996) by analysing the temperature distribution 
of X-ray clusters (Henry & Arnaud 1991). All our models belong to the general class of Cold Dark Matter (CDM) 
scemarios; their linear power-spectrum can be represented as Piin(A;,0) oc k"T^(k), where, for the CDM transfer 
function T(k), we use the Bardeen et al. (1986) flt. In particular, we consider three different EdS models, for which 
the power-spectrum amplitude corresponds to = 0.52 (here ag, is the r.m.s. fluctuation amplitude in a sphere of 
9,h~^ Mpc). They are: a version of the standard CDM (SCDM) model with shape parameter (see its definition in 
Sugiyama 1995) F = 0.45 and spectral index n = 1; the so-called rCDM model, with F = 0.21 and n — 1; a, tilted 
model (TCDM), with n = 0.8 and F — 0.41, corresponding to a high (10 per cent) baryonic content. We also consider 
an open CDM model (OCDM), with matter density parameter Qom ~ 0.3 and erg = 0.87 and a low-density flat 
CDM model (ACDM), with f2om = 0.3, with as — 0.93. Except for SCDM, which is shown as a reference model, 
all these models are also consistent with COBE data; for TCDM consistency is achieved by taking into account the 
possible contribution of gravitational waves to large-angle CMB anisotropics. A summary of the parameters of the 
cosmological models used in this paper is given in Table bJ. 
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Table 1. The parameters of the cosmological models. Column 2: the present matter density parameter Oom; Column 3: the 
present cosmological constant contribution to the density Hqa; Column 4: the primordial spectral index n; Column 5: the 
Hubble parameter h; Column 6: the shape parameter F; Column 7: the spectrum normalization trg; Column 8: the value of the 
parameter 77 in the luminosity-temperature relation required to reproduced the observed log A^-log S (see text for details) . 



Model 






n 


h 


r 




V 


SCDM 


1.0 


0.0 


1.0 


0.50 


0.45 


0.52 


-0.8 


rCDM 


1.0 


0.0 


1.0 


0.50 


0.21 


0.52 


0.0 


TCDM 


1.0 


0.0 


0.8 


0.50 


0.41 


0.52 


-0.3 


OCDM 


0.3 


0.0 


1.0 


0.65 


0.21 


0.87 


-0.3 


ACDM 


0.3 


0.7 


1.0 


0.65 


0.21 


0.93 


-0.2 



4.2 The method 

Theoretical predictions for the spatial two-point correlation function in the RASSl Bright Sample have been here 
obtained in the framework of the above cosmological models by using a method presented in more detail in (Moscardini 
et al. 1999). Here we only give a short description. 

Matarrese et al. (1997; see also Moscardini et al. 1998) developed an algorithm to describe clustering in our past 
light-cone, where the non-linear dynamics of the dark matter distribution and the redshift evolution of the bias factor 
are taken into account. In the present paper we adopt a more refined formula which better accounts for the light-cone 
effects (see Moscardini et al. 1999). The observed spatial correlation function ^obs in a given redshift interval Z is 
given by the exact expression 

Cobs(r) = ^ -2 , (8) 

[f^dziJ^{z,)r{zi)-^] 

where £,obj{r, zi, Z2) is the correlation function of pairs of objects at redshifts zi and 22 with comoving separation r 
and J^{z) is the actual redshift distribution of the catalogue. A related approach to the study of correlations on the 
light-cone hypersurface has been recently presented by Yamamoto & Suto (1999) and Nishioka & Yamamoto (1999) 
within linear theory and by Matsubara, Suto & Szapudi (1997) in the non-linear regime. 
An accurate approximation for ^obj over the scales considered here is 

^ohjir,Zl,Z2) ~ 6cff(zi)fecff (z2)Cm(r, Zavo) , (9) 

where is the dark matter covariance function and Zavc is an intermediate redshift between zi and Z2 , for which an 
excellent approximation is obtained through Z)+(2avc) = D+{ziy^^ D+(z2y^^ (Porciani 1997), with D+{z) the linear 
growth factor of density fluctuations. 

In our treatment we disregard the effect of redshift-space distortions. Some analytical expressions have been 
obtained in the mildly non-linear regime, by using either the Zel'dovich approximation (Fisher & Nusser 1996) or 
higher order perturbation theory (Heavens, Matarrese & Verde 1998). The complicating role of the cosmological 
redshift-space distortions on the evolution of the bias factor has been considered by Suto et al. (1999). A rough 
estimate of the effect of redshift-space distortions can be obtained within linear theory and the distant-observer 
approximation (Kaiser 1987; see Zaroubi & Hoffman 1996 for an extension of this formalism to all-sky surveys). In 
this case the enhancement of the redshift-space averaged power spectrum is given by the factor 1 + 2/3/3 + /3^/5, 
where /3 ~ Q,Q'^/bcg and feoff is the effective bias (see below). Plionis & Kolokotronis (1998), by analysing the XBACs 
catalogue and using linear perturbation theory to relate the X-ray cluster dipole to the Local Group peculiar velocity, 
found (3 ~ 0.24±0.05. Adopting this approach, Borgani, Plionis & Kolokotronis (1999) conclude that the overall effect 
of redshift-space distortions is a small change of the correlation function, which expressed in terms of ro corresponds 
to an ~ 8 per cent increase. 

The effective bias fecff appearing in the previous equation can be expressed as a weighted average of the 'monochro- 
matic' bias factor b{M, z) of objects of some given intrinsic property M (like mass, luminosity, ...), as follows 

h,«{z)=U{zy^ dlnM' b{M',z) Af{z,M') , (10) 



where J^{z, M) is the number of objects actually present in the catalogue with redshift in the range z, z + dz and 
M in the range M, M + dM, whose integral over In Af is N{z). In our analysis of cluster correlations we will use for 
N{z) in eq.(^) the observed one, while in the theoretical calculation of the effective bias we will take the N{z, M) 
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predicted by the model described below. This phenomenological approach is self-consistent, in that our theoretical 
model for N{z, M) will be required to reproduce the observed cluster abundance and their log A'^-log S relation. 

For the cluster population it is extremely reasonable to assume that structures on a given mass scale are formed 
by the hierarchical merging of smaller mass units; for this reason we can consider clusters as being fully characterized 
at each redshift by the mass M of their hosting dark matter haloes. In this way their comoving mass function n{z, M) 
can be computed using an approach derived from the Press-Schechter technique. Moreover, it is possible to adopt for 
the monochromatic bias 6(M, z) the expression which holds for virialized dark matter haloes (e.g. Mo & White 1996; 
Catelan et al. 1998). Recently, a number of authors have shown that the Press-Schechter (1974) relation does not 
provide an accurate description of the halo abundance both in the large and small-mass tails (e.g. Sheth & Tormen 
1999). Also, the simple Mo & White (1996) bias formula has been shown not to correctly reproduce the correlation of 
low mass haloes in numerical simulations. Several alternative fits have been recently proposed (Jing 1998; Porciani, 
Catelan & Lacey 1999; Sheth & Tormen 1999; Jing 1999). In this paper we adopt the relations recently introduced by 
Sheth & Tormen (1999), which have been shown to produce an accurate fit of the distribution of the halo populations 
in the GIF simulations (Kauffmann et al. 1999). They read 



n(z,M) 
and 

b(M,z) -- 



2a 3HSO.on 



TT 8-kG MD+{z)aM 



D+{z)aM 



2p' 



dlnuM 



dlnM 



exp 



2Dl{z)al 



aS'i 



iDliz) 



l + [V^S,/(aMD+{zWP 



(11) 



(12) 



Here o-fj is the mass- variance on scale Af, linearly extrapolated to the present time {z = 0), and Sc the critical linear 
overdensity for spherical collapse. Following Sheth & Tormen (1999), we adopt their best-fit parameters a = 0.707, 
p = 0.3 and A ~ 0.3222, while the standard (Press & Schechter and Mo & White) relations are recovered for a — 1, 
p = and A = 1/2. Notice that 

-1/2 



J\f{z,M) =4Tvr''{z)^ 
dz 



l + ^non r\z) 



n{z,M) <t>(z,M) 



(13) 



where ^(z, M) is the isotropic catalogue selection function, which also accounts for the catalogue sky coverage, as 
detailed below. 

The last ingredient entering in our computation of the correlation function is the redshift evolution of the dark 
matter covariance function ^m- As in Matarrese et al. (1997) and Moscardini et al. (1998) we use an accurate method, 
based on the Hamilton et al. (1991) ansatz, to evolve into the fully non-linear regime. In particular, we use the 
fitting formula given by Peacock & Dodds (1996). 

In order to predict the abundance and clustering of X-ray clusters in the RASSl Bright Sample we need to 
relate X-ray cluster fluxes into a corresponding halo mass at each redshift. The given band flux S corresponds to 
an X-ray luminosity Lx = 4,-Kd\S in the same band, where dL = {1 + z)r{z) is the luminosity distance. To convert 
Lx into the total luminosity Lboi we perform band and bolometric corrections by means of a Raymond-Smith code, 
where an overall ICM metallicity of 0.3 times solar is assumed. We translate the cluster bolometric luminosity into a 
temperature, adopting the empirical relation T = A I/boi (1 + ^)~^ where the temperature is expressed in keV and 
I/boi is in units of 10** /i"'^ erg s~^. In the following analysis we assume A = 4.2 and B = 1/3; these values allow a 
good representation of the local data for temperatures larger than « 1 keV (e.g. David et al. 1993; White, Jones & 
Forman 1997; Markevitch 1998). Analysing a catalogue of local compact groups, Ponman et al. (1996) showed that at 
lower temperatures the Lboi — T relation has a steeper slope [B ~ 0.1). For these reasons we prefer to fix a minimum 
value for the temperature at T = 1 keV. Moreover, even if observational data are consistent with no evolution in the 
Lhoi — T relation out to z ~ 0.4 (Mushotzky & Scharf 1997), a redshift evolution described by the parameter rj has 
been introduced to reproduce the observed log A'^-log S relation (Rosati et al. 1998; De Grandi et al. 1999) in the 
range 2 x 10"^* < S < 2 x 10"^^ (see also Kitayama & Suto 1997; Borgani et al. 1999). The values of rj required 
for SCDM, rCDM, TCDM, OCDM and ACDM models are reported in Table |. A general discussion of the effects 
of different choices of the parameters entering in this method (e.g. the scatter in the Lboi — T relation) is presented 
elsewhere (Moscardini et al. 1999). 

Finally, with the standard assumption of virial isothermal gas distribution and spherical collapse, it is possible 
to convert the cluster temperature into the mass of the hosting dark matter halo, namely (e.g. Eke, Cole & Frenk 
1996) 



T ■- 



7.75 / M 

/3tm I IQisft-iMe 



2/3 



(1 + ^) 



^0 



Q.^{z) 



1/3 



Avi,(2) 

178 



1/3 



(14) 
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The quantity Avir represents the mean density of the virialized halo in units of the critical density at that redshift 
(e.g. Bryan & Norman 1998 for fitting formulas). We assume /3tm = 1.17, which is in agreement with the results 
of different hydrodynamical simulations (Bryan & Norman 1998; Gheller, Pantano & Moscardini 1998; Frenk et al. 
1999). 

Once the relation between observed flux and halo mass at each redshift is established we account for the 
RASSl Bright Sample sky coverage r2sky(>S') (see Figure 7 in De Grandi et al. 1999) by simply setting 4n(l}{z,M) = 

aky[S'(z,M)]. 

4.3 Results 

In Figure ^ we compare our predictions for the RASSl spatial correlation function in different cosmological models to 
the observational data. All the EdS models here considered predict too small an amplitude. Their correlation lengths 
are smaller than the observational results: we find ro ~ 11.5, 12.8, 14.8 Mpc for SCDM, TCDM and tCDM, 
respectively. On the contrary, both the OCDM and ACDM models are in much better agreement with the data and 
their predictions are always inside the 1-cr errorbars (ro — 18.4, 18.6 Mpc, respectively). In order to quantify the 
differences between the model predictions and the observational data, we use again the maximum likelihood approach. 
The minimum value for S is obtained for the ACDM model. A similar value, corresponding to As = 0.7, is obtained 
for the OCDM model, while for the EdS models the resuffing S are much larger: As = 10.6, 19.3, 26.8 for tCDM, 
TCDM and SCDM models, respectively. 

To evaluate what is the effect of neglecting the description of clustering in the past light-cone, as usually done 
in previous analyses, we estimate the cluster correlation function as ^(r) = b'^ff{zincd)^m{r, Zmcd), where Zmcd is the 
median redshift of the catalogue. For the RASSl Bright Sample we have Zmcd — 0.08. The resulting values of the 
correlation lengths obtained in this way are ro ~ 13.5, 15.5, 18.2, 21.9, 22.4 /i"^ Mpc, for SCDM, TCDM, tCDM, 
OCDM and ACDM, respectively. They are typically 20 per cent higher than the estimates obtained by our method. 
This difference is due to the fact that in the past light-cone formalism the matter correlation functions (and bias 
factors) are weighted by a factor N{z) /r{z), for which the average value on the whole sample corresponds neither to 
the median nor to the mode of the redshift distribution. Actually, the presence of the comoving radial distance r{z) 
at the denominator tends to shift the "effective redshift" to smaller values of z. As a consequence, the true cluster 
correlations, which are indeed measured in our past light-cone, have typically smaller amplitudes than those estimated 
at the median redshift of the catalogue. Of course, the importance of this effect becomes larger when deeper surveys 
are considered. 

In order to study the possible dependence of the clustering properties of the X-ray clusters on the observational 
characteristics defining the survey, we use our model to predict the values of the correlation length ro in catalogues 
where we vary the limiting X-ray flux Sum or luminosity Lum (both defined in the energy band 0.5 - 2 keV). Notice 
that this analysis can be related to the study of the richness dependence of the cluster correlation function. In fact, a 
change in the observational limits implies a change in the expected mean intercluster separation d^. Bahcall & West 
(1992) found that the Abell clusters data are consistent with a linear relation ro = 0.4dc, while a milder dependence 
is resulting from the analysis of the APM clusters (Croft et al. 1997). Our results are shown in Figure |^. All the 
cosmological models display a similar trend: in the flux and luminosity intervals here considered, the correlation 
length ro has a slow growth with S'nm (left panel), and a more marked one with Lum (right panel). For example, for 
OCDM and ACDM the correlation length changes from ro — 18 to ro — 21h~^ Mpc, when the limiting flux varies 
from Slim = 10~^^ to Su^n = 10"^^ erg s ^ cm , and from ro ~ 18 to ro — 30/i ^ Mpc, when the limiting luminosity 
varies from Lum = lO**^ to Lum ~ lO***/!"^ erg s~^. The values of ro for the EdS models have similar variations but 
are always smaller. We can compare these predictions to the results obtained by computing the two-point correlation 
function in the RASSl Bright Sample with the same cuts in X-ray flux or luminosity. The estimates of ro obtained 
in this way are presented in Table ^ (where we also report the number of clusters inside each subsample) and shown 
in the figure (open squares). In both panels the whole catalogue is represented by the square on the right. Because of 
the small number of clusters in these subsamples, we prefer to fit the correlation function by fixing the value of the 
slope 7 = 1.8; we find that our values of ro are only slightly dependent on this assumption. The errorbars shown in 
the figure correspond to an increase As = 4 with respect to the minimum value of S. With the assumption that As 
is distributed as a distribution with a single degree of freedom, this corresponds to 95.4 per cent confidence level. 
By analysing the trend with changing limiting flux, we find that the observed values of ro are almost constant even if 
5iim changes by a factor larger than 2. This result is in agreement with what Borgani, Plionis & Kolokotronis (1999) 
obtained for XBACs. We notice that OCDM and ACDM are able to reproduce the amount of clustering shown by 
the RASSl Bright Sample, while all the EdS models strongly underpredict the amplitude of the correlation function. 
The situation is slightly different when we analyse luminosity-limited catalogues. The RASSl Bright Sample suggests 
a small increase of ro with Lum, even if the hypothesis of a constant correlation length cannot be rejected. This is 
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10 100 

r (Mpc/h) 

Figure 3. Comparison of the observed spatial correlation for clusters in the RASSl Bright Sample (shown as in Figure |2|) with 
the predictions of the various theoretical models: SCDM (solid line), rCDM (dotted line), TCDM (short-dashed line), OCDM 
(long-dashed line) and ACDM (dotted-dashed line). 



Table 2. The correlation length ro in flux-limited (left) and luminosity-limited (right) subsamples of the RASSl Bright Sample. 
Column 1: the limiting X-ray flux S'um (in units of 10""'^'^ erg s~^ cm~^) or luminosity Liim (in units of 10'**/i~^ erg s~^); 
Column 2: the number of clusters in the subsample; Column 3: the correlation length ro (in units of Mpc) as obtained 
from the maximum likelihood analysis with fixed slope 7 = 1.8 (errorbars correspond to 95.4 per cent confident level) 



•Slim 


no. of clusters 


'■0 




no. of clusters 


ro 
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27.2t^;5 
30.2l^i 
32. 5l^-^ 
32.8l^J 



consistent with a similar analysis made by Abadi, Lambas & Muriel (1998) on the XBACs catalogue. As shown in 
the left panel of Figure ^, our models always tend to predict smaller correlations, even if the non-EdS models are still 
marginally consistent with the observational data. 



5 CONCLUSIONS 

In this paper we have studied the two-point correlation function ^(r) of a flux- limited sample of X-ray galaxy clusters, 
the RASSl Bright Sample. These observational results have been used to test a theoretical model predicting the 
clustering properties of X-ray clusters in flux-limited surveys in the framework of different cosmological scenarios. 
Our main results are: 
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log Sji^ (erg s-1 cm-2) log L^^ (h-^ erg s'^ 



Figure 4. The behaviour of the correlation length ro as a function of the limiting X-ray flux Su^-^ (left panel) and luminosity 
^lim (right panel). The open squares and errorbars (95.4 per cent confidence level) refer to different subsamples of RASSl 
Bright Sample here analysed. The prediction of the various theoretical models are shown as in Figure^: SCDM (solid line), 
tCDM (dotted line), TCDM (short-dashed line), OCDM (long-dashed line) and ACDM (dotted-dashed line). 



• Assuming an Einstein-de Sitter model, ^(r) can be well fitted using the standard power-law relation ^ = {r/ro) 
with ro = 21.5tlpi-~^ ^pc and 7 = 2.11+° j^ (95.4 per cent confidence levels with one fitting parameter). The values 
obtained in models with matter density parameter Qom ~ 0.3 are quite similar. 

• The amplitude of the correlation function is almost constant when the RASSl Bright Sample is analysed with 
different limiting fluxes in the range 5*11111 = 3 — 7 x 10"^'^ erg s~^ cm~'^ while it displays a slightly increasing trend 
when computed in catalogues with an increasing X-ray luminosity limit in the range Lum ~ 0.01 — 0.4 x 10*''/i~^ erg 
s-\ 

• The comparison with the predictions of our theoretical models shows that the Einstein-de Sitter models are 
unable to reproduce the observational results obtained for the whole RASSl Bright Sample. On the contrary, good 
agreement is found for the models with matter density parameter Qom ~ 0.3, both with and without a cosmological 
constant. 

• Our models are also able to reproduce the behaviour of ro with Sum and Lum, but the observed amount of 
clustering is reproduced only by the open and A models. 

In conclusion, we believe that the method presented here leads to robust predictions on the clustering of X-ray 
clusters; its future application to new and deeper catalogues will allow to provide useful constraints on the cosmological 
parameters. 
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